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Understanding E2 versus 5^2 Competition under Acidic 
and Basic Conditions 

Lando P. Wolters,''^ Yi Ren and F. Matthias Bickelhaupt*^"''^ 



Our purpose is to understand the mechanism through which 
pH affects the competition between base-induced elimination 
and substitution. To this end, we have quantum chemically in- 
vestigated the competition between elimination and substitu- 
tion pathways in H2O + C2H5OH2+ and OH" + C2H50H, that is, 
two related model systems that represent, in a generic 
manner, the same reaction under acidic and basic conditions, 
respectively. We find that substitution is favored in the acidic 



case while elimination prevails under basic conditions. Activa- 
tion-strain analyses of the reaction profiles reveal that the 
switch in preferred reactivity from substitution to elimination, 
if one goes from acidic to basic catalysis, is related to (1)the 
higher basicity of the deprotonated base, and (2) the change 
in character of the substrate's LUMO from C"^— H bonding in 
C2H5OH2+ to C^-H antibonding in C2H5OH. 



Introduction 



Bimolecular base-induced 1,2-elimination (E2) is one of the 
most elementary reactions in organic chemistry.^^"^^ Typically, in 
E2 reactions, an anionic base abstracts a proton from the 
p-carbon center of a substrate molecule while, simultaneously, 
a leaving group at the a position is released, as shown in 
Scheme 1 a. Elimination reactions are, in principle, in direct 
competition with another textbook organic reaction, namely, 
bimolecular nucleophilic substitution (Sn2).^^^ In the case of the 
reaction shown in Scheme 1 a, the competing Sn2 reaction re- 
sults from the anionic reactant carrying out a nucleophile 
attack at the a-carbon center, leading to release of the same 
leaving group as in the E2 reaction. E2 elimination has been 
observed to be favored by stronger bases. In particular, acid 
catalysis, which involves a comparatively weak, neutral base or 
nucleophile, generally goes with substitution, whereas basic 
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Scheme 1. Model E2 and Sn2 reactions corresponding to a) basic and 
b) acidic conditions. 



conditions, often featuring stronger, anionic bases, tune reac- 
tivity towards elimination. 

Elimination and substitution reactions have been studied 
also extensively in the gas phase using mass spectrometry 
techniques.^^"^^^ Note that differentiation between E2 and Sn2 
processes of the anionic reaction systems through mass spec- 
trometry is problematic, requiring additional techniques or 
a special design of reactants (e.g., leaving group connected to 
substrate via a second bond) as both pathways lead, in princi- 
ple, to the same detectable anionic product, namely, the leav- 
ing group Y" (see Scheme 1 a). This problem does not occur 
when studying the corresponding reactions after protonation 
of nucleophile/base and leaving group (Scheme 1 b).^^^'^^'^^^ The 
resulting cationic reaction systems lead to different product 
ions for the respective mechanistic pathways, namely the con- 
jugate acid XH2^ for E2 and the substituted substrate for Sn2, 
which in general have a different m/z ratio that allows for 
straightforward characterization of the active mechanisms. 

Theoretical studies on both anionic and cationic E2 and Sn2 
reactions have provided detailed information on reaction po- 
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tential energy surfaces (PES) and the structure of species and 
transition states (TS)J^^"^^^ Activation-strain analyses^^^^ of the 
anionic substitution reactions between halides and halome- 
thanes^^^^ have shown that the Sn2 barrier decreases as the 
highest-occupied molecular orbital (HOMO) of the nucleophile 
becomes a better electron donor, that is, as the basicity of the 
nucleophile becomes stronger, because of a more stabilizing 
interaction between nucleophile and substrate. Likewise, 
a weaker bond between carbon and the leaving group (C— LG) 
also leads to a lower Sn2 barrier because the energetic strain 
associated with such a weaker bond is less destabilizingJ^^'^^^ 

Herein, we wish to gain insight into the electronic mecha- 
nism behind the observed shift from substitution to elimina- 
tion if one goes from acidic to more basic conditions. To this 
end, we have analyzed the competition between antiperipla- 
nar 1,2-elimination and backside substitution in the cationic 
and anionic model systems H20 + CH3CH20H2^ [Equations (1)] 
and OH" + CH3CH20H [Eqs. (2)], using density functional theory 
(DFT). Our model systems represent, in a generic manner, one 
particular reaction under acidic and basic conditions, respec- 
tively. 

H2O + CH3CH2OH2+ ^ H3O+ + CH2=CH2 + H2O (1a) 
H2O + CH3CH2OH2+ ^ CH3CH2OH2+ + H2O (lb) 
OH" + CH3CH2OH ^ H2O + CH2=CH2 + OH- (2a) 
OH" + CH3CH2OH ^CH3CH20H + OH" (2b) 

Our computations show that indeed substitution is favored 
in the acidic case whereas elimination prevails under basic 
conditions. Activation-strain analyses of the reaction profiles 
reveal that the switch in preferred reactivity from substitution 
to elimination, if one goes from acidic to basic catalysis, is di- 
rectly related to the significantly higher proton affinity of the 
anionic base. Deprotonation enhances the attack of the base 
in both pathways, elimination and substitution, but this effect 
is counteracted by a stronger C— LG bond in the neutral sub- 
strate. Interestingly, however, our activation-strain analyses 
show that protophilic attack benefits more from increasing the 
basicity than nucleophilic attack due to a different composition 
of the substrate LUMO under different reaction conditions. 

Theoretical Methods 

Computational details 

All calculations were carried out using the Amsterdam Density 
Functional (ADF) program^^°^ and the Quantum-regions Inter- 
connected by Local Descriptions (QUILD) program.^^^^ The nu- 
merical integration was performed using the procedure devel- 
oped by te Velde et al.^^^^ The molecular orbitals (MOs) were ex- 
panded in a large uncontracted set of Slater-type orbitals 
(STOs): TZ2P (no Gaussian functions are involved). The TZ2P 
basis set^^^^ is of triple-^ quality for all atoms and has been aug- 
mented with two sets of polarization functions, that is, 2p and 
3d on H, and 3d and 4f on C and 0. An auxiliary set of s, p, d, 



f and g STOs was used to fit the molecular density and to rep- 
resent the Coulomb and exchange potentials accurately in 
each self-consistent field (SCF) cycle. All electrons are included 
in the variational treatment (no frozen-core approximation 
used). 

Energies and geometries were calculated at the OLYP level 
of the generalized gradient approximation (GGA), which in- 
volves the optimized exchange (OPTX) functional proposed by 
Handy and coworkers,^^"^^ and the gradient-corrected functional 
of Lee, Yang and Parr for correlation.^^^^ Scalar relativistic effects 
were accounted for using the zeroth-order regular approxima- 
tion (ZORA).^^^^ Previous ab initio benchmark studies show that 
the OLYP functional, combined with the TZ2P basis set, leads 
to the same trends and qualitative features of the potential 
energy surfaces as CCSD(T)/aug-cc-pv(T+d)Z.^^^^ 

Equilibrium structures were obtained by optimizations using 
analytical gradient techniques.^^^^ At each step in the optimiza- 
tion, the estimate of the true Hessian is improved by an updat- 
ing procedure using the difference of current and previous 
gradients in relation to the difference in geometries. Depend- 
ing on the type of calculation, different Hessian update 
schemes are used. For equilibrium geometry optimizations and 
intrinsic reaction coordinate (IRC) calculations with ADF, the 
Broyden-Fletcher-Goldfarb-Shanno (BFGS) scheme^^^^ is used. 
For transition state searches with QUILD, a weighted combina- 
tion of Powell symmetric-Broyden (PSB)^'^°^ and Murtangh-Sar- 
gent (Symmetric Rank-One, SRI)^"^^^ is used, as developed by 
Bofill.^"^^^ Energy minima have been verified through vibrational 
analysis.^"^^^ All minima were found to have zero imaginary fre- 
quencies, whereas all transition states have one. The relative 
energies in the present study refer to the electronic energies 
without zero-point energy correction. 

Activation-strain analyses 

Insight into the energetics is obtained through activation- 
strain analyses. The activation strain model^^^^ is a fragment- 
based approach to understand the bonding energy AE of two 
fragments during a chemical process and to explain it in terms 
of the original fragments. Within this approach, the potential 
energy surface AE{Q is decomposed along the reaction coordi- 
nate ^ (or just at one point for example at the transition 
state, TS, where AE^^AE^^js) ii^to the strain energy AE^^,^■^^{Q 
that is associated with the geometrical deformation of the indi- 
vidual fragments as the process takes place, plus the actual in- 
teraction energy AE■^^^{Q between the deformed fragments 
[Eq. (3)]. 

AE(a=AE3train(a+AFi,,(a (3) 

The interaction energy AEint(S) between the deformed frag- 
ments is further analyzed in the conceptual framework provid- 
ed by the Kohn-Sham molecular orbital method.^"^"^^ We find 
that donor-acceptor interactions between occupied orbitals 
on one fragment with unoccupied orbitals on the other frag- 
ment, including the HO MO- LUMO interactions, play an impor- 
tant role in AE^^^ and in determining reactivity trends. The 
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PyFrag program was used to facilitate the analyses along the 
potential energy surfacesJ"^^^ Contour values used for the plots 
of the LUMOs in Figure 1 are ±0.02, ±0.05, ±0.1, ±0.2 and 
±0.5. 



former [Eq. (4a)] and 21.7kcalmol ^ exothermic in the latter 
case [Eq. (46)].^"'^ 

H2O + CH3CH2OH2+ ^ H3O+ + CH3CH2OH (4a) 

OH" + CH3CH2OH ^ H2O + CH3CH2O- (4b) 



Results and Discussion 

Reactant complexes 

Our Z0RA-0LYP/TZ2P computations show that our model elim- 
ination and substitution reactions proceed via the formation of 
reactant complexes that, both under acidic and basic condi- 
tions, can rearrange into various variants. These complexes are 
shown in Scheme 2, which also provides relative energies and 
key geometry parameters of each stationary point (for full 
structural data and energies, see Table S1 in the Supporting In- 
formation). The most stable encounter complex in all four reac- 
tions features a hydrogen bond between H2O [Eqs. (1 a) and 
(1 b)] or OH" [Eqs. (2 a) and (2 b)] and a hydrogen atom of the 
leaving group of the substrate. An obvious side reaction for 
our model systems R— 0H2^ and R— OH is deprotonation of the 
leaving group, which is 22.0kcalmor^ endothermic in the 
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Scheme 2. Relative energies (in kcalmor^) and selected bond lengths (in A) of isomeric reactant complexes for 
elimination and substitution reactions of a) H2O + C2H5OH2+, and b) OH +C2H5OH, computed at Z0RA-0LYP/TZ2P 
(see also Scheme 3). 



Note however that this side reaction is not of interest in the 
present work, which serves to model the competition between 
elimination and substitution pathways of the more common 
and ubiquitous reactions of ethers R— OR', amides R— NR'R" and, 
in general, substrates R— L with a leaving group that is not 
easily deprotonated under the reaction conditions. In the fol- 
lowing, we therefore focus on the elimination and substitution 
reaction pathways [Eqs. (1) and (2)]. 

In the acidic case [Eq. (1a,b)], the most stable hydrogen- 
bonded reactant complex is bound by -22.8 kcalmol"^ and di- 
rectly connected to the Sn2 transition states (Scheme 2 a). On 
the other hand, along the E2 pathway, this reactant complex 
first undergoes a change in conformation via a low-barrier 
(0.9 kcalmol"^) rotation of the OH-H2O moiety around the C"- 
0 bond to yield a slightly (0.4 kcal mol"^) less stable species. 
The latter again directly leads to the E2 transition state 

(Scheme 2a). Under basic condi- 
tions [Eq. (2a,b)], the lowest- 
energy reactant complex is even 
more stable, featuring a hydro- 
gen-bond strength of —41.2 kcal 
mor\ Prior to entering into the 
elimination and substitution 
channels, the base has to mi- 
grate first from the leaving 
group to a p hydrogen, yielding 
pre-reactive complexes at —15.9 
and —14.5 kcalmor\ respective- 
ly (Scheme 2 b). Essentially, these 
complexes differ only in the ori- 
entation of leaving group, and 
are easily interconverted via rela- 
tively low-barrier rearrange- 
ments. 

Reaction profiles 

The key stationary points and 
transition states (TS) for all four 
reaction paths, as they emerge 
from our Z0RA-0LYP/TZ2P com- 
putations, are schematically dis- 
played in Scheme 3, including 
relative energies and selected 
5,^2 geometric data (for full structural 

and energy data, see Table S2 in 
the Supporting Information). 
Here we note that Scheme 3 
shows the lowest-energy path- 
ways. Alternative pathways exist 
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0.0 -14.5 +9.1 -14.5 0.0 

Scheme 3. Relative energies (in kcalmor^) and selected bond lengths (in A) of stationary points along the lowest energy elimination and substitution path- 
ways of H2O + C2H5OH2+ [Eq. (1 a,b)] and of OH +C2H5OH [Eq. (2a,b)], computed at Z0RA-0LYP/TZ2P (see also Scheme 2). 



in which stationary points (including the TS) adopt different 
conformations in which one or both of the OH or H2O groups 
are rotated. Typically these are only a few tenths of a kcal mol"^ 
higher in energy and will not be discussed here. For acidic con- 
ditions, we find that elimination [Eq. (1 a)] is outperformed by 
substitution [Eq. (1 b)], with reaction barriers of +7.2 versus 
+ 4.1 kcalmor\ respectively. However, going to basic condi- 
tions, the reaction barrier for elimination [Eq. (2 a)] becomes 
—7.0 kcal mol"^ and thus drops below that for substitution 
[Eq. (2 b)], which assumes the significantly higher value of 
+ 9.1 kcal mol"^ (see Scheme 3). Thus, our model systems 
nicely reproduce the experimental observation that the pre- 
ferred pathway often shifts from substitution to elimination if 
one goes from acidic to basic conditions. Before we provide 
detailed activation-strain analyses of this phenomenon, we will 
discuss in the following the mechanistic features and reaction 
profiles in more detail. 

Under acidic conditions, we find that the elimination reac- 
tion proceeds in a concerted manner, via an E2 mechanism 
which is characterized by a double-well potential energy sur- 
face. One TS, at +7.2 kcalmol"^ relative to reactants, separates 
the aforementioned hydrogen-bonded reactant complex at 
—22.4 kcalmol"^ from the product complex at —11.6 kcalmol"^ 
(see Scheme 3). The product complex consists of a H30^-H20 
moiety interacting with ethylene. Dissociation yields the prod- 
uct olefin +H3O+-H2O at -1.1 kcalmol"^ relative to the reac- 
tants, which makes the overall elimination process slightly exo- 
thermic (see Scheme 3). The competing substitution also pro- 
ceeds via a double-well potential energy surface, as is 
common for gas-phase Sn2 substitution at carbon.^"^^"^^^ The 
barrier for substitution is +4.1 kcalmol"^ relative to reactants, 
that is, approximately 3 kcalmol"^ lower than that for the elim- 



ination pathway (see Scheme 3). Overall, the symmetric b^z re- 
action is thermoneutral. 

Under basic conditions, we find again a double-well poten- 
tial energy surface for the symmetric, thermoneutral Sn2 reac- 
tion, with a TS at +9.1 kcalmor\ that is, only slightly higher 
than for the corresponding substitution pathway in the acidic 
case. At variance, the TS of the competing alkaline elimination 
reaction drops significantly and, interestingly, the elimination 
reaction now proceeds via an Elcb mechanism. Thus, we arrive 
at a triple-well^^^^^ potential energy surface, along which a 
P-proton transfer from substrate to base occurs in a separate 
reaction step prior to the expulsion of the leaving group. The 
occurrence of an Elcb mechanism is in line with the strong ba- 
sicity of the hydroxide anion and the poor leaving-group abili- 
ty of the hydroxyl group.^"^^^ The first TS, associated with 
P-proton abstraction, is at —14.2 kcal mol"^ relative to reac- 
tants, that is, less than 2 kcalmol"^ above the preceding reac- 
tant complex. The resulting intermediate, at —15.2 kcal mol"^ 
relative to reactants, is a 0— H-C" hydrogen-bonded complex 
between the newly formed water and the "C'^H2CH20H carban- 
ion. The second and overall TS, at —7.0 kcal mol"^ relative to 
reactants, leads to the expulsion of the leaving group. We note 
that the Elcb potential energy surface around the TS for leav- 
ing-group expulsion is extremely flat and may be considered 
a transition plateau, in line with earlier 2D explorations of po- 
tential energy surfaces of anionic E2 reactions.^^^^^ This causes 
the position of this TS along the reaction coordinate (but not 
its energy) to be sensitive for algorithmic differences (different 
Hessian update schemes, see computational details) between 
TS optimization and intrinsic reaction coordinate (IRC) explora- 
tion. Thus, the TS optimization yields a saddle point at a stretch 
in C-LG bond of slightly less than 1 A (open dot in PES for 
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basic elimination in Figure 1 b), whereas the IRC scan reaches 
its highest point at a stretch of more than 1.5 A (filled dot in 
PES for basic elimination in Figure 1 b). Importantly, the differ- 
ence of less than 1 kcalmol"^ between the two energies (—7.0 
and —6.1 kcalmor\ respectively, relative to reactants) is signifi- 
cantly smaller than the effect that we are studying and there- 
fore does not affect our conclusions. In the course of its expul- 
sion, the leaving group forms a strongly hydrogen-bonded 
H0"-H20 moiety that interacts with ethylene, at —27.7 kcal 
mol"^ relative to reactants. Dissociation yields the product 
olefin + OH -H2O at -21.8 kcal mol"^ relative to the reactants. 
This makes the overall elimination process under basic condi- 
tions rather exothermic. 

In the following section, we use the activation-strain model 
to show how the barriers arise, how their relative order de- 
pends on nucleophilicity and leaving-group ability, and why 
the preferred reactivity shifts from substitution to elimination if 
one goes from acidic to basic conditions. 

Activation strain analyses: Acidic conditions 

The Z0RA-0LYP/TZ2P activation-strain diagrams (ASD) for all 
four model reactions are shown in Figure 1 a,b. Herein, the 
bold lines represent the reaction profile, that is, the relative 
energy profiles AE{Q as a function of the IRC projected onto 
the stretch of the bond between the a-carbon atom and the 
leaving group (C-LG). A dot on these lines indicates the TS. 
The relative energy is decomposed into a strain energy term, 
^^strain®A whlch Is the dashed line above AE=0, and the inter- 
action energy AE^^^iQ, represented by the dashed line below 
AE=0 (see [Eq. (3)] and Theoretical Methods for details). 

Three general trends emerge from the ASDs of our model 
reactions: In the first place, the strain and the interaction 
curves show counteracting trends, if we go from acidic to 
basic conditions. Thus, for both elimination and substitution, 
the strain curve becomes more destabilizing and the interac- 
tion curve more stabilizing. This leads to a partial cancellation, 
and, as a result, changes in the total energy profile are compa- 
ratively small. Secondly, however, the changes in strain and in- 
teraction, if we go from acidic to basic conditions, are signifi- 
cantly larger for the elimination pathway. Thus, while in the 
acidic case the strain and interaction curves do not differ much 
for elimination and substitution, they adopt significantly larger 
absolute values for elimination under basic conditions. And 
third, the single most determining factor working in the direc- 
tion of the overall shift in preferential reactivity from substitu- 
tion in acidic to elimination in basic conditions is the enor- 
mous strengthening in base-substrate interaction in the latter 
case. In the following, we address the electronic mechanisms 
behind these trends. 

The ASD for the acidic case shows that the barrier for the 
elimination reaction is not only higher than for the substitution 
reaction (+7.2 kcal mol"^ versus +4.1 kcal mol"^) but also 
occurs at a much later stage (see Figure 1 a). This is a conse- 
quence of the fact that, although the proton abstraction and 
leaving-group expulsion occur in a concerted manner, they do 
not occur synchronously: leaving-group expulsion is ahead of 



proton abstraction. One factor working in this direction is a rel- 
atively weakly pulling base, H2O (with a proton affinity of only 
164.5 kcal mor\ calculated from enthalpies at 298.15 K and 
1 atm), and an excellent leaving group, also H2O, which makes 
for a rather soft C— LG bond. A deeper reason for the late, E1- 
like elimination mechanism is revealed when we inspect the 
Kohn-Sham orbital electronic structure of the substrate. For 
this purpose, we describe the molecular orbitals (MOs) of the 
CH3CH2L substrate (L= ^0H2) in terms of a methyl and a substi- 
tuted methyl radical, *CH3((3) + *CH2L(a) (see Figure 1 g, a and 
p label the methyl fragment containing C" and C"^, respective- 
ly). Under acidic conditions, the substrate is protonated at the 
leaving group, and this stabilizes the corresponding *CH2L(a) 
fragment molecular orbitals (FMOs) considerably. 

This has an important consequence. The empty e^(a) FMO 
drops so much in energy that it begins to interact predomi- 
nantly with the occupied e(p) instead of its direct counterpart 
e''(p), and the LUMO of the CH3CH2L substrate becomes mainly 
the antibonding e(p)-e^(a) combination (see Figure 1 g). Thus, 
the LUMO keeps its regular C"-L antibonding character stem- 
ming from e^(a),^^^''^ but it becomes C"^- H bonding (instead of 
C^-H antibonding) due to e(p). This can also be nicely seen in 
the quantitative contour plots of the CH3CH2OH2+ LUMO in 
Figure 1 e, in particular, as one approaches the TS. Thus, the in- 
teraction with the base does not favor C^-H bond rupture, 
and p-proton abstraction occurs late when the C"-L stretch is 
advanced and quite some strain has already been built up. Fur- 
thermore, once the C'^-H bond begins to stretch, this deforma- 
tion has an energy-raising effect on the LUMO (because this 
orbital is C^-H bonding) and thus counteracts the energy-low- 
ering effect of C"-L dissociation. This is at variance with the sit- 
uation for the Sn2 pathway along which the LUMO energy 
drops more significantly because only the C"-L has to break. 
The earlier and more pronounced decrease in LUMO energy 
along the Sn2 pathway is plotted in Figure 1 c. All together, 
these effects favor an earlier onset of the base-substrate inter- 
action along the Sn2 than along the E2 pathway (see Fig- 
ure 1 a). The result is the late. El -like E2 elimination with 
a somewhat higher barrier than the corresponding Sn2 substi- 
tution under acidic conditions. 

Activation strain analyses: Basic conditions 

The ASD for basic conditions shows the aforementioned am- 
plification of destabilizing strain and stabilizing interaction, es- 
pecially for elimination, as compared with the ASD for acidic 
conditions (see Figure 1 b). The increase in the strain curve 
along the Sn2 pathway is caused by the fact that the C"— OH 
bond in ethanol is significantly stronger (with a homolytic 
bond dissociation energy, BDE, of +93.7 kcal mol"^) than the 
C"-0H2+ bond (heterolytic BDE= +32.7 kcalmol"'). The stron- 
ger interaction results from the higher-energy HOMO of OH", 
which makes this species a more strongly binding electron-do- 
nating agent than H2O in the interaction with an electron-ac- 
cepting substrate.^^^'"^^^ The increase in interaction would have 
been even higher if the substrate would not also change from 
protonated to non-protonated ethanol. The neutral substrate 
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Figure 1. Activation-strain analyses of model reactions (a, b; dot designates transition state, TS), substrate-LUMO 
energies along reaction coordinate (c,d), substrate-LUMO contour plots (e,f) and schematic substrate-LUMO com 
position (g,h) under acidic and basic conditions, based on Z0RA-0LYP/TZ2P computations (see text). 



in tine basic case lias a liiglier-energy LUMO tlian tine former in 
wliicli all orbitals are lowered in energy due to the net positive 
potential. The two effects nearly cancel and, therefore, the Sn2 
reaction barrier of +9.1 kcalmol"^ under basic conditions is 
not very different from the Sn2 barrier of +4.1 kcalmol"^ 
under acidic conditions. 



Before continuing with the 
ASD for the elimination process, 
we first elaborate on the nature 
of the substrate's LUMO. The 
latter undergoes a qualitative 
change in character if we go 
from the acidic to the basic sit- 
uation, and this change plays 
a key role in understanding the 
shift in reactivity towards elimi- 
nation. Again, we describe the 
MOs of the CH3CH2L' substrate 
(this time, L' = OH) in terms of 
a methyl and a substituted 
methyl radical, •CH3(P) + -CH2L'(a) 
(see Figure 1 h). Under basic con- 
ditions, the substrate is no 
longer protonated at the leaving 
group and, as a result, the corre- 
sponding *CH2L'(a) FMOs are 
considerably higher in energy 
than in the acidic case. Thus, the 
e^(a) approaches e^(|3), its coun- 
terpart on *CH3(p), from below 
and the overall LUMO of the 
substrate CH3CH2L' becomes 
mainly e^(p) + e^(a), that is, the 
bonding combination of the two 
empty methyl FMOs. Therefore, 
under basic conditions, the sub- 
strate LUMO is o* antibonding in 
both the C"-LG and the C^-H 
bond (see Figure 1 h). This can 
also be clearly recognized in the 
quantitative contour plots of this 
LUMO as the substrate adopts 
the geometry that it has in the 
TS for elimination (see Fig- 
ure 1 f). 

Now, we return to the ASD. 
The combined C"— LG as well as 
C"^— H antibonding character of 
the substrate LUMO leads to 
a completely different behavior 
along the elimination pathway 
under basic conditions. In this 
case, the interaction with the ap- 
proaching base OH" induces 
a weakening not only in the C"— 
LG bond but also in the C^-H 
bond, especially after this bond has been slightly extended 
(see "LUMO in TS" in Figure 1 f which also holds for slightly 
stretched C'^-H distances). Furthermore, the stretching of the 
C'^— H bond that goes with the elimination pathway now con- 
tributes to a lowering of the LUMO energy which therefore 
drops much faster than along the Sn2 pathway, that is, exactly 
the other way around than under acidic conditions (cf. Fig- 




•CH2L'(a) 
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ure 1 c,d). On one hand, the simultaneous stretching of C'^— H 
and C"— LG that is now induced by the base-substrate interac- 
tion causes the strain to increase more pronouncedly. Howev- 
er, it also makes the interaction curve gain in stabilization sig- 
nificantly more steeply. Eventually, it is the better electron-do- 
nating capability of OH" (with a proton affinity of 400.7 kcal 
rr\o\~\ calculated from enthalpies at 298.15 K and 1 atm) in 
combination with the rapidly dropping LUMO energy along 
the elimination pathway that makes the base-substrate inter- 
action pull the barrier of this pathway below that for Sn2 sub- 
stitution (cf. Figure 1 a,b). 

Finally, the strong increase in interaction and, in particular, 
the tf-H antibonding character of the substrate LUMO also 
causes p-proton transfer to run ahead of leaving-group expul- 
sion. This is so, even to the extent where p-proton transfer be- 
comes a separate first step in the reaction mechanism, fol- 
lowed by leaving-group expulsion from the carbanion inter- 
mediate in a second step. Thus, our model reactions show 
a shift from an E1-like E2 elimination that is dominated by Sn2 
substitution under acidic conditions, to an E1cb elimination 
that dominates Sn2 substitution under basic conditions. 

Conclusions 

Our computations confirm the often observed shift from sub- 
stitution to elimination when changing from acidic to more 
basic conditions. This follows from our theoretical analyses 
based on relativistic density functional theory (DFT) of the 
mutual competition between base-induced 1,2-elimination and 
nucleophilic substitution in model systems that represent, in 
a generic manner, one particular reaction under the two differ- 
ent conditions, namely, H20 + CH3CH20H2^ for acidic and 
OH" + CH3CH20H for basic conditions. In particular, we find 
that the elimination pathway in our model systems goes from 
an E1-like E2 mechanism that is dominated by S^2 substitution 
to an E1cb mechanism that prevails over Sn2 substitution. 

Our activation-strain analyses (ASA) reveal that the dominant 
cause for the above shift from nucleophilic (substitution) to 
protophilic reactivity (elimination) is an enormous gain in stabi- 
lizing interaction between the reactants if one goes to basic 
conditions, in particular in the case of elimination. The en- 
hanced interaction is a direct consequence of the fact that the 
base (or nucleophile) changes from H2O to OH". The latter has 
a significantly higher-energy HOMO and thus enters into more 
stabilizing interactions with a substrate LUMO. In the case of 
Sn2 substitution, the enhanced interaction is approximately 
canceled by the more destabilizing strain that arises as we go 
from a relatively weak carbon— leaving-group bond under 
acidic (C"— 0H2^) to a stronger one under basic conditions 
(C"-OH). Thus, the S^2 barrier changes comparatively little. 

On the other hand, the enhancement of the interaction, 
from H2O to OH", is significantly more pronounced for elimina- 
tion and pulls the corresponding barrier below that for substi- 
tution. Interestingly, there is a fundamental reason why elimi- 
nation benefits more from a higher-energy HOMO of the base 
under basic conditions, namely, the fact that also the substrate 
LUMO changes its character such that it favors the elimination 



pathway when it interacts with an attacking base. Thus, the 
LUMO goes from C"^— H bonding and C"— LG antibonding under 
acidic conditions (i.e., in CH3CH20H2^) to antibonding in 
both bonds, C"^— H and C"— LG, under basic conditions (i.e., in 
CH3CH2OH). Therefore, under acidic conditions, the HOMO- 
LUMO interaction between base and substrate does not assist 
tf-H rupture and the LUMO is not stabilized as the tf-H bond 
eventually breaks. At variance, under basic conditions, the 
HOMO-LUMO interaction between base and substrate does 
assist tf-H rupture and the LUMO is stabilized as the C^^-H 
bond begins to stretch, right from the onset of the elimination 
process. 
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